Spectral Densities of Response Functions for the 0(3) Symmetric Anderson and Two 

Channel Kondo Models 



S.C. Bradley 1 , R. Bulla 2 , A.C. Hewson 1 , and G-M. Zhang 3 
department of Mathematics, Imperial College, London SW7 2BZ, UK. 
2 Theoretische Physik III, Elektronische Korrelationen und Magnetismus, Universtitat Augsburg, 86135 Augsburg, 

Germany 

3 Centre for Advanced Study, Tsinghua University, Beijing 100084, P.R. of China 

The 0(3) symmetric Anderson model is an example of a system which has a stable low energy 
marginal Fermi liquid fixed point for a certain choice of parameters. It is also exactly equivalent, 
in the large U limit, to a localized model which describes the spin degrees of freedom of the linear 
dispersion two channel Kondo model. We first use an argument based on conformal field theory to 
establish this precise equivalence with the two channel model. We then use the numerical renormal- 
ization group (NRG) approach to calculate both one-electron and two-electron response functions 
for a range of values of the interaction strength U. We compare the behaviours about the marginal 
Fermi liquid and Fermi liquid fixed points and interpret the results in terms of a renormalized Ma- 
jorana fermion picture of the elementary excitations. In the marginal Fermi liquid case the spectral 
densities of all the Majorana fermion modes display a |cj| dependence on the lowest energy scale, 
and in addition the zero Majorana mode has a delta function contribution. The weight of this delta 
function is studied as a function of the interaction U and is found to decrease exponentially with U 
for large U. Using the equivalence with the two channel Kondo model in the large U limit, we deduce 
the dynamical spin susceptibility of the two channel Kondo model over the full frequency range. We 
use renormalized perturbation theory to interpret the results and to calculate the coefficient of the 
lnu; divergence found in the low frequency behaviour of the T — dynamic susceptibility. 

I. INTRODUCTION 

The 0(3) symmetric Anderson model is a modified form of the symmetric Anderson model which was introduced 
by Coleman and Schofielca. The Hamiltonian can be expressed in the form, 
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+ U{n dr -^)(n dl -±), (2) 

where c da and c da are creation and annihilation operators for the localized impurity d state with spin cr, with Vj = V, 
hj = (Vq + V)/2 and V2 = (Vq — V)/2. The electrons in the impurity state interact via a Coulomb matrix element 
U , and are coupled by the hybridization terms to the conduction electrons which are in the form of a tight-binding 
chain with c\ l(j: c na , the creation and annihilation operators at site n, and t n is the nearest neighbour hopping 
matrix element. The model is basically a symmetric Anderson model (V — Vq) with an anomalous hybridization 
term (V ^ Vo) which breaks both spin and charge conservation. It takes a simpler form when expressed in terms of 
Majorana fermions operators. For the impurity d electrons the Majorana fermion operators are defined by 



d 3 = ~^( c d4 + c di) d o = ^( c dl ~ C dl)' ( 4 ) 
which satisfy the commutation relations, 

{d a ,d/s} = S a ^, (5) 
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where {} indicates an anticommutator. 

The Majorana operators for the conduction electrons are similarly defined by 



X1 („) = _((-1)X T +4 T ), 

■ in-rt/2 

X2 (n) = ^-((-l)»c nT -4 T ), 

_pinn/2 

X3(n) = -^-((-irc w ; + 4;) 

■in-rr/2 

Xo(n) = — 7 =-((-irc nl -cl l ), (6) 

with commutation relations as in (|^). 

Using the Majorana fermion representation the 0(3) Anderson model can be expressed in the form 

3 

H = iJ2 V Xa (0)d a + tV QX o(0)do 
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+ i ^2 ^tnXa{n + l)Xa{n) + Ud 1 d 2 d 3 d , (7) 

a =0n=0 

which is clearly invariant under the 0(3) group of transformations of the basis of the 1,2,3 Majorana fermions. 

There are two aspects of this model that make it worthy of study. One is that the model displays noruFermi 
liquid behaviour corresponding to marginal Fermi liquid theory in lowest order perturbation theory for Vb = 03. The 
marginal behaviour results from singular -scattering with an uncoupled local Majorana fermion mode. Numerical 
renormalization group (NRG) calculationsH have established that this marginal Fermi liquid behaviour persists down 
to T = and that there is an associated residual entropy of ln2/2. This fixed point behaviour has been interpreted 
in terms of the free renormalized Majorana fermion modes. The many-body excitations constructed from these free 
Majorana fermions account for the excitations found at the fixed point in the NRG calculations! For Vq ^ there is 
no uncoupled local Majorana fermion mode and the thermodynamic behaviour corresponds to Fermi liquid theory in 
terms of renormalized Majorana fermions. n 

The original reason Coleman and Schofieldtl studied this model is because, in the large U regime under a Schrieffer- 
Wolff transformation, it maps into a localized spin model, the a-T modelcl, which has a similar form to the two channel 
Kondo model (TCKM). The difference is that the localized spin S — 1/2 is coupled via an exchange interaction to 
the spin and the isospin of conduction electrons in a single channel rather than the spinsj-pf the conduction electrons 
in two distinct channels as in the TCKM (for a comprehensive review of this model seeQ). The Hamiltonian of the 
model can be written in the form, 

oo 

H = S d ■ [.h<j(Q) + J 2 f(0)] + Y, *»(4+i,«r<W + cl a c n+ha ), (8) 

n— 0,(7 

where 5<j is the operator for the localized spin and Ji, J 2 are the two exchange couplings. 
The spin operators, 

*■»(») = \( c U C *t - c ll C nl)> ( 9 ) 

in the basis spanned by the two singly occupied fermion states, |0(t), 1(1)) and |1(T), 0(|)), give a representation for 
the SU(2) algebra for spin a =\- The isospin operators, 

r+(n) - (-l) n 4A> r "W = Mr c n|C„ T , 

Tz(n) = -(c^c n1 +cl lCnl -l), (10) 
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give a representation of the same algebra in the space spanned by the zero and double occupation states, |0,0) and 
|1, 1). By applying the Schrieffer- Wolff transformation it can be shown that the o-r model corresponds to the localized 
or large U limit of the 0(3) Anderson model with the coupling J\ and J2 given by 

2V(V + V ) 2V(V - V ) 
Ji = v J 2 = v • (11) 

The channel isotropic model = J2 = J, corresponds to Vq = in the 0(3) Anderson model. In this case one of the 
Majorana fermions modes is uncoupled and the singular scattering of the conduction electrons gives rise to a marginal 
Fermi liquid fixed point. In terms of Majorana fermions the o-t model (|l^) takes the form, 

3 00 . j 

H = iJ2J2 tn ^ n + ^XaO*) - ■ * (0) X * (0) ' (12) 



a =o 71 = 



where x(n) = (Wn), xa{n)> Xain)). 

Coleman et aH put forward arguments that this one band localized model has the-same fixed point and low energy 
behaviour as the two channel Kondo model. A bosonization approach by SchofieldQ supported this conjecture, and 
so did a comparision of our numerical renormalization group resultsa for the large U 0(3) Anderson model with the 
results for the two channel model. 

In the first section of this paper we show that if wc assume linear dispersion for the conduction electrons the 
relationship between the ct-t model and the spin degrees of freedom of the two channel model becomes an exact one, 
and what is more it is not confined purely to the asymptotic low energy regime but applies to all relevant energy 
scales. In the subsequent sections we extend our earlier numerical renormalization group calculations^ for the 0(3) 
Anderson model to the calculation the dynamical response functions at T — both for the marginal Fermi liquid and 
Fermi liquid cases, and interpret these results using renormalized perturbation theory for the Majorana fermionsHu. 
Finally we use the equivalence of the 0(3) Anderson model and the two channel Kondo model in the large U limit to 
deduce the form of the dynamic spin susceptibility for the two channel Kondo model over all relevant energy scales. 
We use renormalized perturbation theory to deduce the coefficient of the lnw divergence found in the low frequency 
behaviour of the real part of the dynamic susceptibility at zero temperature. 



II. EXACT EQUIVALENCE TO THE TWO CHANNEL KONDO MODEL 

We begin by showing that the Majorana fermion description, used to formulate the 0(3) Anderson model and ct-t 
models, can emerge quite naturally for the two channel model as a representation of the algebra of the total spin 
current of the two channels. For linear dispersion and no cut-off we can express the Hamiltonian of the channel 
isotropic TCKM model in the form, 

H = H + Rj 

r°° dx ■■ 4>) w^-.^) ; 

Z7T j=ltr=U J-oo 

Hi= J « S dJs(0), (13) 

a— x : y ,2 

on taking a Fourier transform to a continuum variable x, where Vf is the Fermi velocity, j is a channel index, 
and the conduction electron operators between the colons have to be to normal ordered. We have assumed s-wave 
scattering only and replaced the incoming and outgoing waves with two left-moving electron fields ipj, a (x); Jg(x) are 
the conduction electron spin current operators 

s a being spin- 1/2 matrices. We can also introduce charge and flavour currents 

3,0 

jf(x) = ]T : i>lMtlri>j'A*) ■■ ■ (15) 
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where t a - ■, are generators of an SU(2) symmetry group. Following Affleck and LudwigS, the free part of the Hamil- 
tonian can be rewritten as a sum of three commuting terms by the usual point-splitting procedure (Sugawara con- 
struction) : 



H = 



2vr 



+00 



dx 



- : J c {x)J c {x) : +^ ■ Jf(x) ■ Jf(x) 
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+- : J s (x) ■ J s (x) 



(16) 



while the interaction term is expressed in terms of the electron spin currents and the impurity spin only. The 
information about the number of channels is contained in the commutation relations obeyed by the spin currents 



[J s Q (x), J b s (x')} = ie abc J«(x)6(x - x') + ^S a . b 6'(x - x') 



(17) 



indicating that Jg(x) form an SU(2) level k = 2 Kac-Moody algebra. Meanwhile, the charge and flavour currents 
satisfy 



[J c (x),J c (x')]=2ki5'(x-x'), 
[J](x), J){x')] = ie abc J?(x)6(x - x') 

+ ^S a , b S'(x- x'). 

They form a U(l) Kac-Moody and another SU(2) level-2 Kac-Moody algebra, separately. 

It is now quite natural to introduce a Majorana representation of the spin current operators in the form, 



(18) 



jy(x) = 



-i : xi(x)x2(x) 



(19) 



where Xi{ x ), X2{x), and X3( x ) are left-moving free Majorana fermion fields. It can be shown that this representation 
reproduces the SU(2) level-2 Kas-Mpody commutation relations. Our approach differs from earlier conformal field 
theory approaches to the TCKMlIinla, where Majorana fermions are introduced only at a later stage after bosonization. 
It is important to note that this Majorana representation of the spin currents is only appropriate for the two channel 
model as it leads to a level-2 algebra. It would be inappropriate for the single channel Kondo model where the 
corresponding spin current generates a level- 1 algebra. 

In a similar way, we can also introduce Majorana representations for the flavour currents 



Jf(x) = ~i ■ X2(x)x' 3 (x) :, 
Jf{x) = -i ■ X 3 (x)xi{x) ■■, 
Jf(x) = ~i ■ x'i(x)X2(x) 

which reproduces the commutation relations satisfied by the flavour currents, and 

J c (x) = -2i : x'a{x)x's{x) ■ 



(20) 



(21) 



can represent the charge current operator. Note that x' a with a = 1,2,3,4,5 are also left-moving free Majorana 
fermion fields. It is well-known that the dynamics of charge, flavour, and spin are completely determined by the 
commutation relations of the current operators. Though the spin currents of the two channel Kondo model can be 
represented in terms of three Majorana fermion fields Xa{x) (a = 1,2,3), we emphasize that they can not be given 
any simple physical interpretation in terms of the original conduction electrons tpj ta (x). 

At this point we have the current operator terms in the Hamiltonian as quartic in the Majorana fields. The Sugawara 
construction enables one to write kinetic energy terms, which are quadratic in field operators, as quartic terms. This 
is what was done earlier in writing the free part of the Hamiltonian in form of eqiiation ( [l6| ) , and it is convenient if 
one is pursuing a purely algebraic approach as used in the conformal field theoryo. However for our purposes it is 
more convenient now to perform an inverse Sugawara construction by the usual point-splitting pmcedure again, and 
rewrite the terms quartic in the Majorana fermions as kinetic energy terms which are quadratidlinlj: 
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J c (x)J c (x) :=4J2 ■ X' a (id x ) x ' a (x); 

3 

J f (x)-J f (x) :=2J2--Xa(idx)x' a (x); 

a=l 
3 

J s (x) ■ J s {x) := 2^ : Xa (id x ) Xa (x). (22) 



a=l 

The model Hamiltonian is transformed and divided into the following two parts, 



^ + = / E / dx: X' a (x){id x ) x ' a (x) :, 



-y^. : x(0) x x(0) : . (23) 

H c +Hf describes the non-interacting charge and flavour degrees of freedom. It has a symmetry group U(1)®SU (2)2 = 
50(5) and is expressed by five free Majorana fermion fields x' a ( x ) ( a — 1; 2, 3, 4, 5). H s is the main part of the model 
and describes the spin degrees of freedom with three left-moving Majorana fermion fields Xa {01 = 1, 2, 3) interacting 
with the impurity spin. It has the symmetry SU (2)2 or 50(3) so that the full Hamiltonian has the symmetry group 
50(5) ® 50(3), which is represented by eight different Majorana fermion fields. 

In the two channel model Hamiltonian, H s given in equation (24) is the only part which includes an interaction 
with the impurity spin. This part of the Hamiltonian is exactly equivalent to the vector part of the a-r model defined 
in equation (|l2|). The zero Majorana fermion contribution to ( |l2| ) is completely decoupled from the other modes so 
we can separate it out, and if we take the continuum limit with linear dispersion equation ( |l2| ) can be written in the 
form, 

H = H sc + H vec 
vt f + °° 

H sc = — / dx : xo{x)(id x )xo{x) ■ 



2n 



—00 

3 f + 



H vcc = 7T- / ] / dx '■ Xa{x)(id x )Xa{x) 



-y5 d -:x(0)xx(0):. (24) 

and can identify H vcc with H s . This implies that the application of the <j-t model is not restricted to the very low 
energy regime but can be used to calculate the impurity contribution to the thermodynamics of the two channel 
Kondo model over the full temperature range. This result is exact subject only to the requirement of linear dispersion 
for the conduction electrons (this is not evident in the approaches that use bosonizationLfEJ~li3) . The result implies 
that the 0(3) Anderson model in the large U regime can be used to calculate the spin correlation functions of the two 
channel Kondo model over the complete temperature and frequency range, u>, T << D, where D is the cut-off imposed 
on the conduction electrons. We exploit this equivalence to calculate the dynamic susceptibility of the TCKM in the 
last section of the paper. Though we have restricted our derivation to the channel isotropic case J\ = J2 = J, the 
result can be generalized to the channel anisotropic case J\ ^ J2. 

The fact that we have been able to establish the exact equivalence the linear dispersion TCKM and tr-r models, 
not just for the low energy regime but over the full parameter range of the model, is at first sight a little surprising. A 
feature of the general TCKM is that there is a strong coupling regime.in which the impurity can be overscreened, and 
this strong coupling fixed point was shown by Nozieres and Blandinll3 to be an unstable one. Impurity overscreening 
on the other hand cannot occur in the strong coupling regime of the o~-t model as the spin and isospin are mutually 
exclusive channels and they cannot both screen the impurity simultaneously. The strong coupling limit for the cr-r 
model, in apparent contrast to that of the two channel model, gives a stable fixed point. However, to access the 
overscreened states of the standard TCKM requires the coupling J to be much greater than the bandwidth D. In 
the linear dispersion two channel model D — > 00 so these states are not accessible in this version of the two channel 
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model. We must distinguish two stong coupling limits. One in which the limit J — > oo is taken with a finite cut-off 
or before the limit D — > oo is taken, and one in which the limits are taken in the opposite order. The first leads to 
overscreening and an unstable fixed point whereas the latter leads to a screened stable fixed point. 

A parallel situation arises in the standard Anderson model. If one assumes linear dispersion then the parameter 
regime in which the impurity level td lies outside the conduction band cannot be accessed. The model with the level 
within the conduction band, and the one with the level below the conduction band, in the Anderson model case have 
the same fixed point but they do have different expressions for the Kondo temperature (seea). 

Before leaving this topic we briefly show that, because the Majorana representation of the TCKM is equivalent to 
the large U limit of the 0(3) Anderson model, we can use results for the 0(3) model to confirm thepd^edption of 
the TCKM fixed point as one corresponding to renocmalized asymptotically free Majorana fermionsOO'u'liSlij. Our 
earlier numerical renormalization group calculationsS established that the fixed point of the 0(3) Anderson model 
corresponds to free Majorana fermions independent of U. This implies that the fixed point corresponds to the non- 
interacting model U = 0. The many-body excitations and their degeneracies at the fixed point, as found in the 
NRG calculations, can be explained in terms of the many body excitations built up from the three non-interacting 
hybridized 1,2,3 Majorana fermion modes, combined with the excitations arising from the unhybridized zero Majorana 
fermion modes described by H sc . This interpretation of the fixed point in terms of the 0(3) Anderson model is simpler 
than the equivalent interpretation as the strong coupling limit of the er-r model. In the latter approach a change of 
boundary condition for the 1,2,3 Majorana fermion modes, which interact with the impurity, has to be invoked as a 
result of the strong coupling, relative to that of the uncoupled Majorana fermion mode. This is not necessary in 
the Anderson model description as there is naturally a difference of boundary condition at impurity in the hybridized 
and unhybridized terms. The leading (dangerously) irrelevent interaction term, which leads to the logarithmically 
divergent contributions to the specific heat coefficient and susceptibility in these models, can be interpreted as the 
renormalized interaction term U in the Anderson model (seeati for further details). 

To derive a complete description of the many-body excitations and their degeneracies at the fixed point for the 
TCKM from the U = 0(3) Anderson model we simply have to replace the uncoupled Majorana degree of freedom 
described by H sc by the five uncoupled Majorana degrees of freedom in H c and H{, which describe the charge and 
flavour excitations. There are two slightly different sets of excitations, one set corresponding to a chain of even sites 
(sector A) and another set corresponding to an odd number of sites (sector B). When these are combined, as shown in 
Table 1, they reproduce the many-body energy levels and degeneracies found in the NRG calculations!!-!! and conformal 
field theory calculations!^ for the TCKM. Since the interaction is restricted to the spin part of the Hamiltonian H s 
the leading (dangerously) irrelevent interaction which leads to logarithmically divergent contributions to the specific 
heat coefficient and susceptibility in the TCKM is exactly the same as the renormalized interaction U ~ Ik of the 
0(3) Anderson model, where Xk is the Kondo temperature. We will use this explicitly later in section 4 to calculate 
the dynamic susceptibility of the TCKM from that of the 0(3) Anderson model. 



In this section we extend the numerical renormalization group method to calculate the dynamic features of the 0(31 
Anderson model. For the basic details of the NRG method as applied to this model we refer to the earlier paperEl 
where the thermodynamic properties were calculated. Calculation of the dynamic behaviour via the NRG approach 
is rather more involved because additional matrix elements have also to be calculated in the sequence of iterative 
diagpnalizations . The general details of the NRG approach to the calculation of spectra for impurity models are given 
inE3'Eil Because both charge and spin are not conserved for this model the calculations here differ from the standard 
case in a number of respects and the relevant details are given in the appendix. 



III. NUMERICAL RENORMALIZATION GROUP CALCULATIONS 



A. Spectral densities of single particle Green's functions 



The spectral density A ct (oj) of the single particle Green's function of spin a at the impurity site is given by 




nm 




(25) 
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where E n are the exact eigenvalues of the Hamiltonian and Z is the partition function (Z = J2 n ex P( — PE n ))- This 
is evaluated using the many-body states calculated in the sequence of iterative diagonalizations of the NRG. Because 
of the decreasing energy scales in the sequence of steps in the NRG, the spectral density for a particular frequency 
have to be calculated from the step in the sequence corresponding to the appropriate energy scale. The spectrum is 
then built up from the complete set of NRG steps. Though the method has been developed for finite temperature 
situations, the procedure is more straightforward and more accurate for the case T — > 0, and we have restricted 
the calculations here to the low temperature limit. Though a large number of states are retained at each step in 
the NRG calculation (of the order of 300-500 not including degeneracies), the spectrum is discrete and to obtain a 
continuous curve a broadening is introduced which decreases logarithmically with the energy scale. The advantages 
of this approach over alternative methods, such a Monte Carlo, are that it can be used at very low temperatures, can 
be used for any value of U, and the results can be checked against known sum rules and identities. The details of the 
low energy spectrum.pan be calculated to a very high resolution. For further details of the approach we refer to the 
papers cited earlierElEil. 

Due to the anomalous hybridization term, the 0(3) Anderson model does not conserve either total spin or total 
charge. However, the combination of total spin plus total isospin T is conserved and the many-body states and 
the eigenstates \m) can be classified using the quantum-numbers T and T z associated with the total of this angular 
momentum and its z-component ^pte that in referenceu the symbol j was used for this quantity, but here we revert 
to the earlier notation T used intjEa). In evaluating ( p5| ) the sum over the degeneracies associated with T z can be 
done explicitly as explained in the appendix. The single particle Green's functions for the impurity can be expressed 
as linear combinations of the Green's functions for the Majorana fermions, G a p(t — t') = —iO(t — t')({d a (t)dp(t')}) , 
for a,/3 = 0,1,2,3 using equations (0) and In the absence of applied field these Green's functions are diagonal 
(oc 5 a ,p)i and due to the 0(3) symmetry those corresponding to a = 1,2,3 are equal. We can deduce the corresponding 
spectral densities, p a (tj), using equation (|25|). The results for these spectra for Vq = 0, the non- Fermi-liquid case, 
are shown in figures 1-4 for a range of values of the interaction U, and corresponding results for Vq ^ are shown in 
figures 8-10. 

In figure 1 the spectral density for the uncoupled a — Majorana fermion (scalar) is shown together with those 
for the a = 1, 2, 3 Majorana fermions (vector) which are hybridized with the conduction electrons (and which are all 
identical). There is an additional delta function contribution to the a — spectral density which is not shown but 
its weight as a function of U is shown in figure 6. What is clear is the cusp- like peak in the spectral density of the 
scalar Majorana fermion at Fermi-level and the fact that the height of this peak is U dependent. The spectral density 
for the vector Majorana fermions also has a cusp-like peak at the Fermi-level but, in contrast to the scalar case, its 
value at the Fermi level is independent of U. In figure 2 the spectral densities of the scalar (a = 0) Majorana fermion 
are shown on their own with an inset showing the linear behaviour in the neighbourhood of the Fermi level. The 
corresponding plots for the vector Majorana fermions are given in figure 3, which also show linear behaviour in the 
neighbourhood of the Fermi-level. r-,,-. 

We can analyze these results using the renormalized perturbation theory developed in earlier workuEl. The per- 
turbation theory is in terms of renormalized propagators for the Majorana fermions, and renormalized parameters, 
A = zA, Ao = ZqAq, and U = vzo^Toi23(0, 0, 0, 0), where z, zo,are the wavefunction renormalization factors derived 
from the non-singular contributions to the self energies, E(w) and £o(w), associated with the a = 1,2,3 and a = 
Majorana fermions respectively, A a (ui) = %V^po(uj) where po{ui) is the density of states of the conduction electrons, 
and roi23 > ^2 j W3 , u>4,~) , is the irreducible four vertex for the Majorana interaction. The retarded Green's function 
for the Majorana fermions then has the form, 

= i — -(rl ; -is )t : > ( 26 ) 

LU + lA a - Eaa(w) - Saa(w) 

where fla\<jj) and Ea^(w) are the regular and singular renormalized self-energies, and we have taken the wide band 
limit with A a (u>) independent of lo. In the renormalized perturbation expansion there arc countcrterms which have to 
be included to satisfy the renormalization conditions and the Majorana fields are rescaled to absorb the wavefunction 
renormalization factors so that the free propagators in the expansion are l/(u> + iA a ) (seeu for further details). 
For the non-Fermi liquid case Vq = the leading order low frequency correction to the renormalized self-energy 
for the a = 1, 2, 3 Majorana fermions comes from the second order diagram for the singular part & s '(lj) shown 
in figure 5(a). The result is 

Re£«( w )= (JL) cjln(t), (27) 

and 



ttA / VA 
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for T = 0. As [/ — > 0, A Q — > A Q , {7 — » {/ and this result goes over to that of the standard second order perturbation 
theory. For the density of states of the a = 1,2,3 Majorana fermions in the neighbourhood of the Fermi level we find, 



7rA 



_ 2 

1_ f Gs) M + °(- 2ln2 H) 



(29) 



for a = 1, 2, 3. We see from this result that the density of states at the Fermi level is 1/ttA, independent of U, which 
is the same result as for the 0(4) Anderson model, and that the leading correction term is linear in These results 
explain the features seen in the plots of the spectral densities of the vector Majorana shown in figures 1 and 3. The 
linear term in the spectral density at the Fermi level is seen clearly and the value at the Fermi-level agrees with the 
theoretical value 1/irA to within 2%. 

The second order contribution to the renormalizcd self energy £o(w) for the a = Majorana fermion is the same as 
the for the 0(4) model as it only involves explicitly the propagators of the renormalized a — 1, 2, 3 Majorana fermions 
(see figure 5(b)). To lowest order the result for the imaginary part is 



--!(£)'(£) < 3 »» 

for T = 0. The leading term for the real part of the self-energy is linear in u but is cancelled by the counterterm 
contribution to comply with the renormalization conditions. When this is inserted into equation fl2q ) then we find for 
the corresponding spectral density 

p (u)=z S{ U ) + 1 ^ T (-^) +0(c 2 ), (31) 



2?rA VttA, 

The delta function contribution at zero frequency reflects the fact that the zero mode associated with the uncoupled 
a = local Majorana fermion exists down to T = 0, and it is this mode that leads to the anomalous entropy of ln(2)/2. 
Its weight in the spectral density is reduced to zq, which we cannot calculate within the renormalized perturbation 
approach. However, we can calculate it to leading order for small U/irA, and the result is zq = 1/(1 + (U/irA) 2 ). For 
large U/irA, the model can be related to the two channel Kondo model and there is only one low energy scale, which 
is the Kondo temperature. We expect, therefore, in this limit that zq ~ Ik/ A which implies that for large U/irA, 
that zq should tend to zero exponentially in the variable U/irA. This behaviour is consistent with the NRG results 
for the delta function weight plotted as a function of U/irA shown in figure 6. The non-delta function contribution 
has a finite limit as oj — > given in renormalized perturbation theory by the the second term in equation (|3l]) . Unlike 
the spectra for the a = 1,2,3 Majorana fermions the value of this contribution at the Fermi level depends on U 
and monotonically increases with increase of U/irA. We see from equation ( |3l| ) that to lowest order in U/irA it is 
proportional to U 2 /(irA) 3 . As the second order diagram in the renormalized perturbation theory gives the exact 
result to order oj 2 for ImSo(w) there should be no further contributions to the leading two terms in equation (p3). In 
the large U limit we deduced from a comparison of the results with those for the two channel Kondo model (seeo) that 
U/irA — > 1. The height of the peak at the Fermi-level compared with that for the vector term in this limit is then 
given by zq/2z (using A = zA), and provides an estimate of the ratio of the wavefunction renormalization factors of 
the scalar and vector Majorana fermions. 

The correction terms in (|3ll ) arise from terms in the imaginary part of the self-energy which behave as w 4 for small 
u>. The second order diagram in the renormalized perturbation theory gives an contribution uj 4 but there are also 
contributions to this order from higher order diagrams. More important, however, are contributions to the imaginary 
part which behave like w 2 |w|, as they result in corrections to ( |3l| ) which behave like |cj|. Such a term is generated from 
the higher order diagram shown in figure 7. The asymptotic contribution to the spectral density from this diagram 
as lu — > is 

' U \ 4 \u 



- 2z ^[^a)W ^ 

There may be other higher order diagrams which contribute terms to this order u> but checking this is difficult. The 
explicit evaluation of the higher order diagrams is in general complicated and there does not seem to be any simple 
way of seeing whether or not a diagram gives a contribution of this order. It is clear from the results shown in figure 2 
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that there is a net contribution of this order as there is clear |w| dependence in the neighbourhood of the Fermi level 
If U is small compared to ttA then the second order diagram which gives an uj 2 correction to 
over the leading |u;| term which is of fourth order. The second order result corresponding to 



if? 



) should dominate 
is 




2A 2 



(33) 



The uj 2 behaviour of po(u>) in the neighbourhood of the Fermi level in the small U regime is evident in figure 4 which 
gives a plot of po(ui) for U/nA = 0.01. 

For Vo ^ the zero Majorana fermion is no longer uncoupled from the conduction electrons and for U — the 
corresponding spectral density has a resonance width of Ao. Having lost the zero mode in this channel the scattering 
the vector Majorana fermions with the scalar Majorana fermion with the interaction U is no longer singular. For 
Ao/A <C 1, the resonance in the scalar channel is very narrow compared with that of the vector channel, and there 
is a low energy range in which the behaviour is very similar to that found in the marginal case. There is a cross-over 
eventually on reducing the energy scale to a Fermi-liquid form of behaviour. This is apparent in the plots of the 
spectra in figures 8-10 for the case A /A = 10~ 4 . In figure 8 the spectral densities of both the scalar and vector 
Majorana fermions are shown for the same set of values for U/irA as in figure 1. Due to the small value which was 
used for the ratio Ao/A the spectra for the vector Majorana fermions shown in figure 9 are quite similar to the 
corresponding curves for the marginal Fermi- liquid case shown in figure 3. The difference is only revealed when we 
compare the insets in the two figures. On the very low energy scale around the Fermi level there is a distinct flattening 
and rounding of the spectrum seen in the inset in figure 9, consistent with an uj 2 behaviour in contrast to the form \uj\ 
seen in the inset in figure 3 (which persisted down to the lowest energy scales for which the spectra were calculated) . 

If we apply second order perturbation theory for finite Ao then, in the low frequency range uj <C A, the imaginary 
part of the self-energy 

has the form, 



i , w \ Ao, A lu 2 

( a:J •-•H 1 ' Ai 



(34) 



In the range uj > Aq this behaves essentially as \uj\, as in the non- Fermi liquid case (|28|). Only in the very low 
frequency range uj -C Ao do we get the influence of the Fermi-liquid fixed point and the uj 2 behaviour. In the low 
frequency limit uj <C Aq, the spectral density of the vector Majorana fermions is then given by 



7rA 



2AA P 



2A 



+ 11 + ^11^ 



U 
irA 



(35) 



for a = 1,2,3. For large values of U we can apply the renormalized perturbation approach and the corresponding 
result for the spectrum of the vector Majorana fermions in the very low frequency limit uj -C Aq is 



p a (oj) 



1 

TTA 




(36) 



The spectra for the scalar Majorana fermion for Ao ^ 0, however, are very different from those for Ao = as can been 
seen by comparing the results given in figure 10 with those in figure 2. There is no zero mode and the peaks at the 
Fermi level coincide. In the very low frequency range uj <C Ao we find for Po{uj) using the renormalized perturbation 
theory the result, 



1 



ttA 



2A 2 



Ao 
A 



U 

7rA 



(37) 



The value of the density of states at the Fermi level for the results shown in figure 10 are within 2% of the Freidel 
sum rule result l/7rAo for all values of U. The uj 2 behaviour in the very low energy region near the Fermi- level is 
clear in the inset in figure 10. 







B. Dynamic spin and charge susceptibilities 



We use the same approach to calculate the two particle Green's functions corresponding to the dynamic spin and 
charge susceptibilities. The corresponding spectral densities are of the form 

2 



8(u-(E n -E m )) 



x ( e -?E m _ e -0E n } (38) 

with Q = g /lib (ti^t — n d,|)/2 for spin and Q = [rid t -\ + n c i,[ — l)/2 for charge. This is similar in form to that for 
the single particle Green's function given in equation (p5f), which corresponds to Q = c da with a difference in sign 
in the factor in the last bracket due to the use of the commutator rather than the anticommutator for the spin and 
charge Green's functions. Using the procedure outlined in the appendix these can be expressed in terms of reduced 
matrix elements and the degeneracy with respect to T z can be explicitly summed over (see equation d5l]))- Both the 
operators for spin and charge transform as irreducible tensors of the form Qj, and the result in each case gives 



nm 

( e -PE m _ -j3E n 



5{uo - (E n - E m )) 



(39) 



where \ n Q mj is the reduced matrix element (see equation (|5l|)). 

The dynamic susceptibilities for the non-interacting model (U — 0) can be evaluated analytically. There are two 
parts to the d contribution to the spin susceptibility Xs( w )- O ne arises from the bubble diagram, figure 11 (a), where 
the propagators involved are those corresponding to a component of the vector Majorana fermion which we denote by 
Xj v (w) (in which we absorb factors of g^)- This contribution is just half of that obtained for the standard symmetric 
Anderson model with the real part given by 

A 2 



A/47T 

(u 2 + 4A 2 ) 



In 



A 2 



4A 
H tan 

LU 



'(I) 



and the imaginary part by 



Im*7M = 



A/27T 



(lu 2 + 4A 2 ) [ w 



A 



In 1 



u 
A 2 



tan 



(40) 



(41) 



The other contribution, which we denote by %JJ s (iLi;), is from bubble diagram, figure 11(b), in which one of the 
propagators corresponds to a component of the vector Majorana fermion while the other propagator is for the scalar 
Majorana fermion. The real part is given by 



2 A/4tt 



(lu 2 + A 2 ~ A 2 ) 2 +4w 2 A 2 



- A 2 



2tu 2 



In 



A 2 



A 2 



+ 



2An 



-tan 



Ao 



+ {similar term with Aq <-> A}, 



(42) 



and the imaginary part by 



ImxJTM = 



wAA /4tt 



A 2 



A 2 



A 2 ) 2 
+ A 2 



In 



ojAq 
+ A 2 ) 



4w 2 A 2 



tan 



uJ 

A~ 



A 2 



similar terms with Aq «-> A 



(43) 
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In the limit Ao — > A the vector-scalar contribution is the same as the vector-vector contribution and we recover the 
results for the dynamic susceptibility of the symmetric Anderson model. In the limit Ao — ► 0, where we have a zero 
mode, the real part of the vector-scalar contribution develops a logarithmic singularity, 



(44) 



and the imaginary part has a discontinuity at the origin, 



. sgn(aj)A 

Im ^>) = ^VW (45) 

We also have some exact reults for U 7^ from the Shiba relation which relates the imaginary part of x a {u>)/u) for the 
impurity to the real part of Xs(^) m the limit u> — ► 0, 

f hnxsju + i6) \ (Re Xs (0)) 2 (AR , 
hm -° { 2™ | = -GiZT- (46) 

This result, however, is only applicable to the standard symmetric Anderson model where A = A. 

We can use the exact results we have for this model to gauge the accuracy of our NRG calculations. In figure 12 
we plot the imaginary part of Xs{<*>) divided by the square of the real part for various values of U. According to the 
Shiba relation these should all pass through the same point at u> = and this is clearly seen in the results figure 12. 
The actual value at this point is about 8% too high, reflecting errors of a few % in the individual susceptibilities. 
The overall accuracy over the full frequency range can be gauged by a comparison of the NRG results with the exact 



results deduced from equations (40), (|4lj), (|42| ) and ( [13|) for [7 = 0. In figure 13 we plot the NRG results for the real 
part of Xs(w) for various values of Ao, and compare them with the analytic results. It can be seen that the two sets 
of results are in good agreement over the complete frequency range. 

To see how the dynamic spin susceptibility for finite U differs from that for the standard symmetric Anderson model 
we have plotted both the imaginary part and real part in figures 14 and 15 for U/irA = 1.5 and various values of Ao. 
As we reduce Ao the imaginary part changes more and more rapidly at to = until finally for Ao = it has a finite 
discontinuous jump at this point as can be in figure 16, where the results for Ao = are shown for various values 
of U. The peak in real part of the spin susceptibility narrows considerably on reducing Ao as can be seen in figure 
15. To be able to compare the results for various values of Ao the values shown in figure 15 have been multiplied 
by the factor l/x s (0), so it is not so apparent that the narrowing is accompanied by a sharp increase in the absolute 
value at low frequencies. In the limit Ao = the real part of Xs(w) has a singularity at uj = 0. A discontinuity in 
the imaginary part of Xs( w ) f° r = at w = with a logarithmic singularity in the real part is predicted from 
equations ( |45"| ) and ([3]) for U = 0, so it is not surprising that these features persist for small values of U. For large 
U we can make predictions from our results about the behaviour of the spin susceptibility of the corresponding two 
channel Kondo model over the full frequency range. We look at this in detail in the next section. 



IV. SPIN DYNAMICS OF THE TWO CHANNEL KONDO MODEL 



In the large U regime we can apply the mapping the 0(3) Anderson model in the localized limit and the NRG results 
of previous section to deduce the spin dynamics of the two channel Kondo model over the complete frequency range. 
As the mapping is for the spin degrees of freedom of the TCKM, which is the only part with an interaction term, it 
cannot be used to calculate all the relevant properties of the TCKM, only those where the other degrees of freedom 
of the conduction electrons do not contribute. It can be used for example to deduce the the impurity contribution to 
the thermodyamics, such as susceptibility and specific heat, and also for the spin dynamics of the impurity, where the 
uncoupled degrees of freedom factorize and cancel out. This is a distinct advantage for NRG calculations to use the 
equivalent in one channel model because a larger percentage of the states can be retained after each diagonalization 
in the iterative procedure. The number of states that can retained be at each stage is proportional to 1/4 c , where 
N c is the number of channels, so this can be serious limitation when the number of channels is increased. However, 
quite apart from these practical considerations, the relation between the TCKM and-the the 0(3) Anderson model 
provides a simple way of understanding in terms of renormalized Majorana fermionsaa. 

We saw in the previous section that the real part of the dynamic spin susceptibility (Ao = 0) at T = has a 
logarithmic singularity in u> for U = and that in the NRG results this singular behaviour persists for finite K_Such 
a term in lnu has also been found in direct NRG calculations of the dynamic susceptibility for the TCKMEi To 
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analyse the situation in the large U limit to use the mapping to the TCKM it is useful to work with a conserved 
quantity, which in this case is not the spin but T the spin plus isospin. The corresponding local dynamic susceptibility 
we denote by Xt(w), where t = ad + t<j. Because the charge fluctuations of the impurity are suppressed for large U 
this susceptibility becomes identical with the spin susceptibility in the large U limit. We can use the NRG approach 
described in the last section to calculate Xt(u) and verify this explicitly. We find for U/nA = 4 at the most a 5% 
difference between Xt(w) and Xs(^) over the complete low frequency range lo <C U. 

For the large U regime we can apply our renormalizcd perturbation approach to calculate the uj dependence of xt , 
following the same steps we used in our earlier work to calculate its temperature dependence in the zero frequency 
limits. We first calculate the contribution due to the non-interacting renormalized Majorana fermions. This is given 
by the vector- vector bubble shown in figure 11 (a), together with an overall factor of 4 and A replaced by A so that 
for T = 0, 



RextH 



A 



tt(lu 2 + 4A 2 ) 



hi 



A 2 



4A , /u\] 

- tan (a)) 



(47) 



Because the single particle Green's function and the quasiparticle Green's function in the renormalized perturbation 
approach differ by a wavefunction renormalization factor z, one might have expected an extra z 1 multiplying (f47j ) to 
give the two particle Green's function corresponding to the dynamic susceptibility. This is not the case because there 
is a compensating factor associated with the coupling to the external field. The net result is that these two terms 
cancel so that in the quasiparticle Hamiltonian in the presence of a external field the quasi-particles are coupled to the 
field with the same coupling as the bare particles, i.e. there is no g-factor renormalization. In the zero frequency limit 
( [47| ) gives Xt(0) = 1/ttA, which corresponds to zero temperature impurity contribution to the total susceptibility we 
found in our earlier work in the same limit. The leading order correction term to this result arises from a second order 
term in U, shown in figure 17, where the scalar-vector bubble shown in figure 11, gives a lnw contribution similar to 
the InT dependence in the same susceptibility at ui = seeta, so that asymptotically as u> — > we find 



ln(w/A) 

7rA 



(48) 



The correspondence of the calculation of the lnw terms with the calculation of-,the InT terms for the zero frequency 
susceptibility is so precise that we can appeal to arguments in our earlier works to show that higher order lnw terms 
cancel in the same way as the higher order InT terms. 

Though we cannot explicitly carry out the full renormalization programme to calculate the renormalized parameters 
U and A in terms of the bare parameters U and A, we can use the results U/ttA — 1, 7rA = Tk from our earlier 
work, which were obtained by equating the expressions for the temperature dependence of the susceptibility in the 
large U limit. We then find for the asymptotic behaviour of the dynamic susceptibility of the two channel Kondo 
model xkm(w) as to — ► 0, u> = seeQ, so that asymptotically as ui — > we find 



RexKMM = =- (1 - In(w7r/T K )) . 



(49) 



In figure 18 we plot the real part of both xt(w) and Xs(k>) versus ln(w) for U/nA = 4 and the asymptotic form given 
in equation ( p9[ ) is clearly seen in the low frequency limit. We have also calculated Tk both from the limiting form of 
Xt(T) as T ->■ and xtM as u ->■ 0, for U/nA = 4 and V = 0.00141 and find T K = 1.3 ± 0.2 and T K = 1.28 ± 0.1, 
respectively, which confirms the result (H9) to within the error limits of the calculations. 



V. CONCLUSIONS 



We have shown that 0(3) Anderson model in the large U limit can describe precisely the interacting spin degrees 
of freedom of the linear dispersion two channel Kondo model. This extends earlier work in that it establishes that 
this relationship applies over the full energy range of the two channel model and not just in the immediate region 
of the low energy fixed point. It enables us to calculate the impurity contribution to the thermodynamics, and also 
the impurity spin dynamics, for the TCKM from the 0(3) Anderson model in the strong correlation regime. This 
has a number of advantages. The 0(3) model is most conveniently expressed in terms of Majorana fermions and 
the low energy fixed point of this model corresponds to free Majorana fermions whatever the value of the on-site 
interaction U, so the nature of the fixed point is the same in weak and strong coupling. This is similar to the standard 
Anderson model, which has a Fermi liquid fixed point, independent of the value of U. The main difference between 
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the 0(3) and the standard Anderson model, is that 0(3) model has a zero mode in the case that corresponds to the 
channel isotropic TCKM which results in singular scattering of the renormalized Majorana fermions giving rise to 
logarithmically divergent contributions to the impurity specific heat coefficient, susceptibility, and dynamic spin-spin 
response function. This singular low energy behaviour corresponds the marginal Fermi liquid theory that .was put 
forward to describe the anomalous behaviour observed in the normal state of the cuprate superconductorscJ. The 
description of the TCKM in terms of Majorana fermions not only gives a simple interpretation of the many-body 
excitations at the low energy fixed point, both for the channel isotropic and anisotropic cases, but can also be applied 
on all energy scales. 

To find the renormalizations of the Majorana fermions explicitly we use the NRG approach to calculate the Majorana 
spectral densities. We have calculated weight of the zero Majorana fermion mode and have shown that it decreases 
exponentially with U in the large U regime. We have also demonstrated the |w| behaviour of the other Majorana 
fermion modes in the marginal Fermi liquid case as lo — > 0. We have used the same approach to calculate the two 
particle response functions which we have determined over the full frequency range for T = 0. We have exploited the 
mapping between the models to deduce the local dynamic spin susceptibility of the TCKM. Renormalized perturbation 
theory has been used to interpret the results and to calculate the coefficient of the lnw term in asymptotic form of 
the real part of the dynamic spin susceptibility of the TCKM in the limit u> — > 0. 
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VI. APPENDIX 



Here we outline the steps involved in the evaluation of the spectral density given in equation ( p5| ) . The many-body 
states \n) are classified using the quantum numbers T and T z of the total angular momentum (spin plus isospin) 
and expressed in the form \T n , T z<n , w n ), where w n labels any further degeneracies of the state \n). The degeneracies 
associated with the z-component of angular momentum in (]2^) can be summed over explicitly. The creation and 
annihilation operators, c d a and c\ a , first of all have to be expressed in terms of the irreducible tensor operators V* , 
k=l,0, q=l,0,-l. From equation (32) in referenced, 



1 n 



V 1 - r f 
v l - c d. 



(50) 



Then the Wigner-Eckart Theorem, 



(T'X,w'\V Q k \T,T z ,w 



1 



V2T' + 1 
x (T,T z ,k,q\T'X,) 



(T',w'\\V k \\T, W ) 



(51) 



is used to express all the matrix elements in terms of the reduced matrix elements (T", w'\ \V£\ |T, w) which are 
independent of the T z quantum numbers, where (T, T z , k, q\T', T' zl ) are the Clebsch-Gordon coefficients. For example 
for the spin down matrix elements we find 



l 

72 



T^vn i ^z^m i Win 



13 



1 



2(2T„ + 1) 
T n ,W n 



Vo 1 



-Mi 7 



2 



(52) 



For T n ^ T m , the Vg'-term does not contribute and for this term we find 



E 



'41 



2(2T„ + 1) 



Vo 1 



E 



T T 1 



T T 



(53) 



For the case T n = T rn = T we find 



E 



T, w r 



T, W r , 



1 



T, w„ 



(54) 



The reduced matrix elements (T„, w n 



V- 9 
v 



T m , w m ) are calculated from those of the previous iteration in two steps. 

We start from the expression for the reduced matrix element with respect to the eigenstates of the impurity-conduction 
electron chain with N + 1 sites, 



T n , w r , 



V q 

1 



Tm ) Wn 



N+l 



m , i m, 7 Pn 



V q 
1 o 



Tm, > Tm j Pra j 

N+l 



(55) 



in terms of the basis states 



T n ,r n ,p n ) used for this calculation (these are constructed from the eigenstates of the 

' N+l 



Tn,r n ) and the states for the additional site labelled by p n , sec reference!] for further details). The 

' N 



N site chain 

coefficients UT n (w n ,r n ,p n ) are the components of the eigenstate 

matrix elements ^T„,r„,p 
step as follows: 
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/ N+l 

T m , r mi pm ) have to be related to the reduced matrix elements of the previous 
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The matrix elements between all states at intermediate steps have to be calculated. The initial values are 
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The spectral density of the spin up Green's function can be evaluated in a similar way and is found to depend only 
on the matrix element of the tensor operator Vq and not on V^ . This Green's function is the same as that of the 
a = 1 or a = 2 Majorana fermion Green's functions, which are equal and equivalent to the Green's function for the 
a = 3 Majorana fermion as a result of the 0(3) symmetry of the model. We can deduce that Vq is associated with 
the a = 1,2,3 Majorana fermions and V ° with the a = Majorana fermion. As spin down Green's function is the 
half the sum of the Green's functions for the Majorana fermions for a = 0, 3 we can deduce that the term involving 
Vq in the final formula for the spectral density of the spin down Green's function is the spectral density of the a = 
Majorana fermion (apart from a factor 1/2) and, similarly, the term in Vq 1 is the spectral density of the a = 1,2,3 
Majorana fermions. 



TABLE I. Many-body excitation energies and their corresponding degeneracies (dg) of the two channel Kondo model corre- 
sponding to the single particle spectrum at the low-energy fixed point. The energy levels without primes correspond to sector 
A, and those with primes to sector B. 
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El + El/2 
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FIG. 1. The spectral density p (uj) of the Majorana fermion impurity Green's function in the marginal Fermi liquid situation: 
V = 0.01414, V = 0.0 (A = 10" 4 7T, A = 0.0) . 
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FIG. 2. The spectral density po (ui) of the a = Majorana fermion Green's function in the marginal Fermi liquid situation: 
V = 0.01414, V = 0.0 (A = 10~ 4 tt, A = 0.0) . The inset shows the |w| dependence near the Fermi level. 
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FIG. 3. The spectral density pi,2,3 (w) of the a = 1,2,3 Majorana fermion Green's function in the marginal Fermi liquid 
situation: V = 0.01414, Vo = 0.0 (A = 10~ 4 7r, A = 0.0). The inset shows the |w| dependence near the Fermi level, pi, 2, 3 (0) 
is within 2% of the value given by the Friedel sum rule. 
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FIG. 4. The spectral density po (u>) of the a = Majorana fermion Green's function in the marginal Fermi liquid situation: 
V = 0.01414, Vo = 0.0 (A = 10~ 4 tt, A = 0.0), indicating the influence of w 2 corrections for small U/nA. 
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FIG. 8. The spectral density p (to) of the Majorana fermion impurity Green's function in the Fermi liquid situation: 
V = 0.01414, V = 0.00014 (A = IO^tt, A = 10" 8 7r) . 




FIG. 9. The spectral density pi, 2, 3 (w) of the a = 1, 2, 3 Majorana fermion Green's function in the Fermi liquid situation: 
V = 0.01414, Vo = 0.00014 (A = 10~ 4 7r, A = 10~ 8 7r) . The inset shows the region of the spectrum near the Fermi level where 
flattening consistent with to 2 behaviour is observed, pi, 2, 3 (0) is within 2% of the value given by the Friedel sum rule. 
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FIG. 10. The spectral density po (w) of the a — Majorana fermion Green's function in the Fermi liquid situation: 
V = 0.01414, Vo = 0.00014 (A = 10 _4 tt, A = 10~ 8 tt). The inset shows the oj 2 dependence near the Fermi level. p (0) is 
within 2% of the value given by the Friedel sum rule. 
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FIG. 12. (Imx, (w) /w) /2tt (Rex s (0)) 2 for the standard model: V = V = 0.01414 (A = A = 10~ 4 7r). 
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FIG. 13. The real part of the dynamic spin susceptibility \a (w) for the non-interacting model: V = 0.01414 ( A = 10 4 7r) , 
U/ttA — 0.0. The circles denote analytic values. 
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FIG. 14. The imaginary part of the dynamic spin susceptibility \s (w) hi the Fermi liquid situation: 
V = 0.01414 (A = 10" 4 tt) , U/tyA = 1.5. 
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FIG. 15. The real part of the dynamic spin susceptibility \s (w) Ixs (0) in the Fermi liquid situation: 
V = 0.01414 (A = 10" 4 tt) , U/ivA = 1.5. 
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FIG. 16. The imaginary part of the dynamic spin susceptibility \s (w) in the marginal Fermi liquid situation: 
V = 0.01414, V = 0.0 (A = 10" 4 7T, A = 0.0) . 
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FIG. 18. The real part of the dynamic spin susceptibilities xt (w) and \s (w)in the marginal Fermi liquid situation: 
V = 0.01414, Vo = 0.0 (A = 10~ 4 7r, Ao = 0.0) , U/ttA = 4.0, displaying linear dependence on ln(a>) in the low energy re- 
gion. 
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